Rotating quantum turbulence in superfluid 4 He in the T = limit 
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Observations of quantum turbulence in pure superfluid 4 He in a rotating container are reported. 
New techniques of large-scale forcing (rotational oscillations of the cubic container) and detecting 
(monitoring ion transport along the axis of rotation) turbulence were implemented. Near the axial 
walls, with increasing forcing the vortex tangle grows without an observable threshold. This tangle 
gradually develops into bulk turbulence at a characteristic amplitude of forcing that depends on 
forcing frequency and rotation rate. At higher amplitudes, the total vortex line length increases 
rapidly. Resonances of inertial waves are observed in both laminar and turbulent bulk states. On 
such resonances, the turbulence appears at smaller amplitudes of forcing. 

PACS numbers: 67.80.bd, 65.40.-b 



In superfluid 4 He driven to the state of quantum turbu- 
lence (QT), quantized vortex lines (each with circulation 
K = 1.00 x 10~ 3 cm 2 /s) make a dynamic tangle At 
lengthscales smaller than the mean distance between the 
lines £ = L -1 / 2 , where L is their total length per unit vol- 
ume, QT differs markedly from classical turbulence. In 
the zero-temperature limit T = viscosity vanishes, al- 
though finite dissipation [2| is still maintained probably 
due to short-scale deformations of vortex lines (Kelvin 
waves) - the cascade of energy to these lengthscales being 
maintained by line reconnections. One thus might need 
very small amplitudes of forcing to sustain steady tur- 
bulence. As the rate and efficiency of reconnections are 
expected to decrease with the alignment of neighbouring 
lines, to investigate the intrinsic mechanisms of dissipa- 
tion of QT it is instructive to study the case of vortex 
tangles polarized by an imposed rotation. For instance, 
in classical fluids rotation modifies the dynamics of the 
energy cascade 3] and generally opposes turbulence. 

In rotating superfluid, the equilibrium state is that 
with an array of parallel vortex lines of a uniform den- 
sity Lq = 2Qq/k. Low-frequency (w < 2J7o) excitations 
are expected to be inertial waves [4j (in essence, oscilla- 
tions of compression and bending of the vortex array), 
i.e. the same as in classical fluids [1, @. Swanson et al. 
[fj studied rotating superfluid 4 He at high temperatures, 
at which one can destabilize individual vortex lines via 
an axial counterflow of the normal and superfluid compo- 
nents. They observed two different turbulent states that 
appear after exceeding the critical values of the coun- 
terflow velocity. Using numerical simulations and disre- 
garding vortex pinning at the container walls, Tsubota et 
al. interpreted the first transition as the proliferation 
of Kelvin waves on isolated rectilinear vortex lines (high 
degree of polarization, L rs Lq) and the second one as 
the onset of a vortex tangle of low polarization (L ^> Lq) 
caused by reconnections of neighbouring vortex lines. 

The aim of this work was to observe and investigate 
QT in rotating superfluid 4 Hc in the T = limit, driven 



at classical (3> I) lengthscales; for instance, whether the 
critical amplitude of forcing is finite and if there could be 
more than one turbulent regime? The experiments were 
performed in isotopically-pure liquid 4 He at temperature 
T < 0.2 K and pressure 0.1 bar. A cube-shaped volume 
of side d = 45 mm was confined by six square electrodes 
(see side view in inset in Fig. [I] details are in jsj). The 
whole cryostat could rotate around its vertical axis at a 
computer-controlled angular velocity Q(t). 

To force turbulence, an AC component of frequency 
uj and small amplitude Af2 was added to the DC an- 
gular velocity fio of rotation of the cryostat, f2(i) = 
f2o + Af2sign(sin(o;i)). The following mechanisms of tur- 
bulence generation are expected: large-scale vortices are 
created after each stroke due to the flow separation near 
the corners of lateral walls; quantized vortices are agi- 
tated near the axial walls due to surface friction (pinning- 
unpinning); inertial waves induce large-scale AC flow in 
the container (in resonance conditions, its velocity can 
be much greater than ~ dAtt). 

To detect turbulence in the middle of the container, 
we monitored the connectivity of vortex lines in the ax- 
ial direction through measurements of the transport of 
ions trapped on them. Negative ions (electrons in a bub- 
ble state) were injected from the field-emission tip 1 mm 
below the centre of the grid in the bottom plate - by 
applying voltage to the tip Vti P = 525 V (thus inject- 
ing the current of ~ 100 pA that was mainly terminated 
at the injector grid) for 0.1s, and then reverting Vti p 
back to zero. At the opposite wall, the currents sepa- 
rately collected by the collector, J co ii(i), grid together 
with its frame, I gY id(t), and top plate, -/piate(i), were con- 
verted into voltage and recorded. The grid in front of the 
top collector, made of a square mesh of wire of diameter 
w = 0.02 mm with period s = 0.5 mm, was stretched on a 
metal ring-frame of 13 mm i. d. and 19 mm o. d. and was 
separated from the surrounding top plate by a 1 mm gap 
- thus making the radius of its gridded part R\ — 6.5 mm 
and effective radius dividing the grid and outer top-plate 
i?2 ~ 10 mm. The time constant of the recording elec- 
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Time after start of injection (s) 

FIG. 1. (color online) Records of current to the top collector 
Ico\\(t) after a O.ls-long injection, beginning at t — 0, at the 
bottom of the cell - for different amplitudes of agitation AQ. 
Q = 1.5rad/s, u = 1.5rad/s, T = 0.2 K, E = 40V/cm. 
Inset: sketch of a stationary cell with the vortex array (left); 
sketch of vortex lines in the axial region of the cell at moderate 
(center) and high (right) level of forcing. 
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FIG. 2. (color online) Analysis of the records of currents 
to the top collector (Fig.[T]), grid and surrounding plate for 
Qo = lu = 1.5rad/s. Left axis: grid transparency 9 (solid red 
circles), fraction of charge arriving outside the grid a (open 
red squares). Right axis: width At, of collector current pulses 
(A). The vertical lines indicates the critical amplitudes of 
forcing Afi c i « 0.01 rad/s and AS1 C 2 ~ 0.03 rad/s. 



tronics was 30 ms. All charges, currents and voltages are 
quoted with the opposite sign (i. e. assuming the elec- 
tron's charge to be positive). The voltage between the 
top collector and its grid (and top plate), side plates, 
and bottom plate were 10 V, 100 V, and 190 V, respec- 
tively - making the mean axial field in the drift space 
E = 40V/cm. In Fig. Q] we show I co \i(t), all at the same 
fio = 1.5 rad/s but for different amplitudes of forcing AO 
at frequency uj = 1.5 rad/s. With increasing forcing, the 
amplitude I* of the peak of /coii(t), firstly increases and 
then decreases. The peak's width at half-maximum, At, 
only begins to increase at high amplitudes of forcing. 

At T < 0.2 K all negative ions in 4 He are bound to 
vortex lines. During the injection they create a compact 
charged vortex tangle near the tip that spreads through 
the bottom grid. With little forcing, there is an array of 
rectilinear vortex lines everywhere, except for this vortex 
tangle at the bottom that feeds ions into the array. The 
ions' time of flight along straight vortices to the top, at 
the terminal velocity of ~ lOm/s |9|, is just ~ 4 ms. The 
position of the maximum at 0.20 s and width at half- 
height At = 0.21s of I co ii(i) (Fig. [[]) reflect the slower 
dynamics of the evolution of the charged vortex tangle 
near the injector and not the ion emission duration 0.10 s: 
first ions cross the injector grid some 0.05 s after the be- 
ginning of the injection, the maximal current is fed into 
rectilinear lines after further 0.15 s, followed by a nearly- 
exponential tail with the lifetime of some 0.25 s. 

The total charges arriving at different electrodes, Q C oii, 
Qgrid, Qpiate, were obtained by integrating the corre- 
sponding currents. These allowed to quantify the frac- 
tion of ions arriving outside the radius of R2 ~ 10 mm 
("lateral spread"), a = Q p i a tc/(Qcoii + Qgrid + Qplatc) j 
and transparency of the grid assembly, 9 = Q C oii/(Qcoii + 



Qgrid)- The lateral and axial spread parameters a and 
At are indicative of bulk turbulence (at least when they 
grow above their no-forcing values, see below), while the 
grid transparency 9 is a local measure of the state of vor- 
tex lines near the collector grid. When rectilinear vortex 
lines are not perturbed and their density L is compara- 
ble to that of grid cells, s~ 2 , many lines terminate not at 
the collector but at the grid; hence, only a fraction of ions 
reach the collector. This fraction increases with increas- 
ing agitation of vortex lines that can cause intermittent 
reconections of lines between pinned (terminated on grid) 
and unpinned (transparent for ions) states. When these 
reconnections result in a vortex tangle of density L 3> s~ 2 
near the grid, most ions that arrive within the gridded 
radius do make it to the collector, except for those head- 
ing to a grid's wire head-on within the catching diameter 
w' = w + j£, where 7 ~ 1. Thus, if ions arrive uniformly 
within a radius R2 ~ 10 mm or greater, 



(1) 



evaluated for 7 = 1 and L — t~ 2 = 10 4 cm -2 . 

The graphs of 9, At and a versus AO, corresponding 
to the records from Fig.[TJ are shown in Fig.[2J There 
are two limiting regimes and their cross-over between 
the 'critical' amplitudes AO c i and AO C 2 . The first (bulk 
laminar) regime, at < AO < AO c i, in which At and a 
stay unchanged and only 9 grows linearly with increas- 
ing AO, corresponds to local vortex tangles near the top 
and bottom grid co-existing with an array of rectilinear 
vortex lines elsewhere in the bulk. With increasing forc- 
ing the density of the local tangle eventually becomes 
much larger than Lq, after which 9 is expected to satu- 
rate below the geometrical transparency of the grid as- 
sembly 0.64. Above Af2 c i, we do observe levelling-off at 
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FIG. 3. (color online) The amplitude of the collector current 
peak 7* versus AQ for different Qo and uj. The vertical lines 
indicate the critical amplitudes AS7 c i and AS7 C 2. The scale of 
the top horizontal axis (for Qo/uj — 3) is different from that 
of the bottom axis (for Qo/uj = 1) by the factor of 1.73. The 
solid lines guide the eye through the linear part. The values 
of 7*(A57) for Qo = ui — 1.5 rad/s, calculated using Eq.[2]and 
data from Fig. [5] are shown by asterisks. 



9 w 0.26, which, according to Eq.[TJ corresponds to the 
local density L ~ 10 4 cm -2 (while L Q = 3 x 10 3 cm -2 ). 

In the second regime, at AS1 > ASl c2 , most ions ar- 
rive at the collector much later than the peak injection 
at 0.20 s, with a large spread of arrival times At and 
with a larger lateral spread a than for the injection into 
an array of rectilinear vortex lines. This regime cor- 
responds to a vortex tangle filling the whole container 
(bulk turbulence). The cross-over between the local and 
bulk turbulence takes place gradually between ASl cl and 
Af2 C 2. The first critical amplitude, ASl cl , is set by the 
comparison of local I near the collector grid with the 
grid period s. The second one AJl c2 , corresponds to At 
and a becoming clearly greater than Ato and ctq set by 
the injection's intensity and rotation rate SIq. For the 
shown example of Slo — 1-5 rad/s and uj = 1.5 rad/s, 
AO c x ~ 0.01 rad/s and AS7 C 2 « 0.03 rad/s; however, for 
f?o = 1.5 rad/s and ui — 0.21 rad/s, ASl cl w 0.04 rad/s 
and Af2 C 2 ~ 0.12 rad/s. For all studied frequencies uj 
and Oo, we found that AS1 C 2 ~ 3ASl c i. The saturated 
value of 9 was always within 0.25 ± 0.1. With the total 
injected charge Q to t = Qcoll + Qgrid + Qpiatc = 23 pC 
and rotation rate firj = 1-5 rad/s, the low-forcing values 
were At = 0.21s, Co = 0.29. With decreasing Jlo, 6>o 
decreased while Ato and oo increased. 

To investigate the observed regimes at different Slo and 
uj, we monitored the amplitude of the collector peak 7*, 
which is related (asterisks in Fig. [3]) to the quantities 9, 
a and At through 

7* ^ 0A6Q tot 9(l - a)At~\ (2) 

where 0.46 is an empirical factor. At small Ail < AS! cl , 
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FIG. 4. (color online) Top panel: the ratio of the gradi- 
ent go = (7* (AO,) - 7o)/Af2 of the linear part of 7*(Af2) 
to A/j^ ax (r2n) versus the reduced forcing frequency uj/2Qo 
for three rates of rotation fin- go was measured at am- 
plitude AS! = 0.0075 rad/s (except for open blue squares, 
for which AQ = 0.015 rad/s). A/^ (0.5 rad/s) = 1.95 pA, 
A/max(l-0 rad/s) = 5.0 pA, A7* ax (1.5 rad/s) = 6.2 pA. The 
right axis shows the corresponding critical amplitudes AS7 c i 
from Eq.[31 Resonances are indicated by vertical lines. Bot- 
tom panel: g = (I*(AQ) — Iq)/AQ versus uj/2Qo for Qo = 
1.50 rad/s and four different values of Afi. The right axis 
shows the amplitudes AS7 c i = 6.2pA/(? for Qo =1.5 rad/s. 
The horizontal lines indicate the values of Afi c i(<?) = AQ 
beyond which the linear growth I*(AQ) (i.e. g = go) breaks 
down; from top to bottom, AQ: 0.015, 0.030, 0.045 rad/s (the 
colors of lines and symbols for the same AQ are the same). 



7*(Af2) grows linearly I* = 7q + go AS1, due to the grow- 
ing 0(AQ). At larger AJ1 > AO c2 , I* (AO) decreases 
chiefly because of the increasing Ai(Afi). As the maxi- 



mum. 



is just above Afi c i, the gradient of the linear 



part of 7*(Ar2), go, is inversly proportional to AS1 C 



(3) 



where A7^ ax (Slo) = -^max — ' s independent of uj. 

In Fig. [4] (top panel), we plot the dependences on the 
forcing frequency U) of the ratio go/AI^ x for three rates 
of rotation Slo- There are several broad peaks at frequen- 
cies to that are not fixed but proportional to the rotation 
rate Slo- This is substantiated by plotting versus w/2Slo : 
at uj/2flo < 0.6, all three datasets collapse on the unique 
function ASl~ L 1 (w/2Slo)- 

Resonances of inertial waves are indeed expected at 
certain uj/2Q < 1 @. These were calculated for an 
ideal liquid subject to non-slip boundary conditions in a 
rotating cube [10( but generally disagree with the posi- 
tions of our peaks. Whether this disagreement is because 
of different boundary conditions of turbulent superfluid 
in the T = limit or because our cell is not truly cu- 
bic (gaps between square electrodes at the edges of the 
cube and circular holes in the centres of four electrodes 
for injectors and collectors of ions) is an open question. 
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FIG. 5. (color online) Mean values of vortex line density L, 
measured in the transverse direction, versus forcing amplitude 
Af2 for various Qq. Solid lines indicate the no- forcing limit 
I/O = 2Qo/k. The arrows indicate the values of AQ c i and 
AO C 2 = 3Af] c i, where A!l c i = AI* /go for particular values 
of cj/2f2o were calculated using the data for SIq = 1.5 rad/s: 
AI* = 6.2 pA and (see Fig.Hl bottom panel) go = 300pC/rad 
for u/2S!o = 1, go = 600pC/rad for w/2f2o = 0.5 and go = 
300pC/rad for w/2«o = 0.33. 



The non-monotonic behaviour of 7* (An) at Ail > 
Aflci is presented in Fig. 0] (bottom panel). Here we plot 
g = (J* (Afl) - I *)/An. When An > An c i(uV2n ), 
g no longer coincides with the gradient go but falls be- 
low it. For each forcing amplitude An, this limiting 
value of g, go — AI* nllx ({l )/ AQ, is indicated by a hor- 
izontal line; these lines are in a good agreement with 
the observed plateaus ( "chopped-off peaks") in g. When 
An exceeds An c i(w/2no), the peaks in go(u/2ilo) cor- 
respond to troughs in g(u/2Qo). 

We found that not only Afi c i(w, no) stays the same 
if the ratio w/no is kept constant, but all I* (An) for 
the same no collapse on a single curve if plotted ver- 
sus AO/AQ c i(w,0), Fig.El This reflects the fact (Eq.0J) 
that all functions 9(ui, n ), a(w,n ) and Ai(w,n ) col- 
lapse on single curves when plotted this way (this was 
confirmed independently). In the investigated range of 
resonant frequencies, 0.05 < w/2no < 0.5, the peak (on- 
resonance) values of AO c i seem to scale as oc lu~ x / 2 , that 
can be summarized: An c i ~ 0.01(cj/no) -1 / 2 rad/s. At 
troughs, although never farther away from a resonance 
than one linewidth, ASl c i w 0.04 ± 0.01 rad/s is roughly 
independent of the forcing frequency. In the limit of 
high-frequency forcing uj/2Q,q > 1, where inertial waves 
cannot propagate, all An cl keep increasing in the man- 
ner perhaps compatible with the common scaling law 
An cl oc {uj/VLq) 1 / 2 [12], but eventually saturate at differ- 
ent values seemingly proportional to no. At the highest 
frequency studied lu — 10 rad/s (at no = 1.5 rad/s), Afi c i 
was about 10 times greater than that at the strongest in- 
ertial wave resonance (at u/2Hq — 0.49 rad/s). 

To confirm that in a turbulent state the length of 
the vortex tangle L exceeds Lqi we measured its mean 



value along a horizontal axis of the cube via scattering 
of charged vortex rings 11] of average radius R = 1.6 /jm 
sent from the left injector to the right collector. In Fig.[5j 
one can see that at forcing amplitudes smaller than AO C 2 
the vortex line length, measured far from the horizon- 
tal walls, is only slowly increasing from the equilibrium 
L = 2n /K. But at An > An c2 , L(An) L starts to 
grow faster and rapidly becomes much greater than Lq. 
This supports our interpretation that, at An > An c 2, a 
vortex tangle of low polarization develops in the bulk. 

In the limit of slow rotation, the dependence L(An) in 
the turbulent regime (Fig. [5]) is steeper than at An/n <C 
1 and seems to be closer to L oc An 3 / 2 (as for no = 0). 
The exponent 3/2 is indeed expected if we assume that 
the fully-developed turbulence exerts the resistive torque 
oc An 2 ; then its time-averaged work oc An 3 should 
equate the rate of dissipation oc v'(kL) 2 . The critical 
amplitude for sustainable QT at no = is small but 
increases seemingly proportionally to no. 

With fast rotation (no > w/2), inertial waves are an 
inherent feature of the large-scale superfluid dynamics, 
thus greatly reducing the critical amplitude of forcing 
at strongest resonances. Between the resonances, the 
critical amplitudes are roughly frequency- independent, 
An c i w 0.04 rad/s, either due to the dominance of other 
mechanisms of generating turbulence or complex inter- 
action of inertial waves of different wavenumbers. With 
further increasing no, as the on- resonance critical ampli- 
tudes increase as oc (no/w) 1 / 2 , one can expect the iner- 
tial waves to become non-dominant again if the forcing 
frequency oj is kept constant. 

Summary. This is the first experimental studies of 
rotating QT in the T = limit. The new detection 
techniques exploited the quantized nature of vorticity 
and were sensitive to the density and degree of entan- 
glement of vortex lines separately in the bulk and near 
the container walls. Within our resolution, no finite crit- 
ical amplitude of forcing was required to sustain either 
bulk turbulence in a non-rotating container or boundary 
turbulence near the walls of rotating container. In rota- 
tion, with increasing the amplitude of the AC-rotational 
forcing, we observed a gradual growth of vortex tangles 
near the horizontal walls followed by a transformation, 
around the finite critical amplitude of forcing AO C 2, of 
an array of rectilinear vortex lines into a developed bulk 
vortex tangle. Inertial waves, observed at frequencies 
cj/2no between 0.03 and 0.5, help generate turbulence: 
on-resonance, the critical amplitude An c 2(w/2no) can be 
several times smaller than when off-resonance. 

This work was supported by the Engineering and 
Physical Sciences Research Council [grant numbers 
GR/R94855, EP/H04762X and EP/I003738]. 
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